The Crystal Field Potential of PrOs 4 Sbi 2 : Consequences for Superconductivity 
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The results of inelastic neutron scattering provide a solution for the crystal field level scheme in 
PrOs4Sbi2, in which the ground state in the cubic crystal field potential of Th symmetry is a Ti 
singlet. The conduction electron mass enhancement is consistent with inelastic exchange scattering, 
and we propose that inelastic quadrupolar, or aspherical Coulomb, scattering is responsible for 
enhancing the superconducting transition temperature. PrOs4Sbi2 appears to be the first compound 
in which aspherical Coulomb scattering is strong enough to overcome magnetic pair-breaking and 
increase T c . 

PACS numbers: 74.70.Tx,71.27.+a,75.10.Dg,78.70.Nx 



Praseodymium filled skutterudite compounds, with 
general formula PrT4Mi2, where T is one of the tran- 
sition metals Fe, Ru, or Os, and M is a pnictogen (P, As, 
or Sb), show a remarkable variety of interesting physi- 
cal phenomena, including metal- insulator transitions 0, 
quadrupolar heavy fermion behavior 0, 0], and super- 
conductivity Q. In particular, PrOs4Sbi2 has attracted 
attention as a heavy fermion superconductor, in which 
quadrupolar fluctuations may play an important role in 
the pairing mechanism 0. This proposal is based on 
two key observations; (i) there is a significant f-electron- 
induced mass enhancement of the conduction electrons, 
observed in specific heat, upper critical field 0-E1- and 
de Haas- van Alphen measurements , and (ii) the mag- 
netic susceptibility indicates that the crystal field ground 
state is non-magnetic 

A knowledge of the crystal field ground state is essen- 
tial to understanding the role of the f-electrons in the 
superconductivity. In analyzing the magnetic suscepti- 
bility and specific heat, Bauer et al considered two pos- 
sible crystal field models In cubic symmetry, the 
Pr 3+ ion splits into a singlet (Fi), a non-magnetic dou- 
blet ^3), and two magnetic triplets (T4 and T5). Crys- 
tal field models with either the Ti singlet or doublet 
as ground state were both broadly consistent with the 
data; in both cases, the T5 triplet was the lowest ex- 
cited level estimated to be at less than 1 meV in energy. 
A T3 non-Kramers doublet ground state is of particu- 
lar interest as it provides the necessary conditions for 
quadrupolar Kondo fluctuations to be responsible for the 
heavy fermion behavior |8| , and was favored by analyses 
of the entropy 0,^3- However, the alternative Ti sin- 
glet ground state has also been proposed following exper- 
iments that have explored the crossover to a field-induced 
ordered phase 0, 0, ^| , so this important question re- 
mains unresolved. 

Inelastic neutron scattering is the most direct method 



of determining the crystal field potential and level scheme 
of metallic rare earth systems. In this report, we present 
the results of a comprehensive set of measurements of 
crystal field transitions in PrOs4Sbi2 as a function of 
temperature. From a simultaneous profile refinement of 
all the spectra, normalized on an absolute intensity scale, 
we have concluded that the Ti singlet is the ground state 
level. Discrepancies with earlier neutron scattering re- 
ports [3 are explained by the need for extra terms in 
the cubic crystal field Hamiltonian that are required by 
the Th point group symmetry |l4|. Although our data 
are inconsistent with a quadrupolar Kondo scenario, we 
conclude that inelastic quadrupolar fluctuations do play 
a vital role in enhancing the superconducting transition 
temperature, through aspherical Coulomb scattering of 
the conduction electrons |15| . 

We performed our experiments on the same polycrys- 
talline sample that was used for previous neutron scat- 
tering measurements; details of its preparation and char- 
acterization can be found in Ref. |6J. The inelastic neu- 
tron scattering experiments have been performed on the 
time-of-flight Fermi chopper spectrometer LRMECS at 
the pulsed spallation neutron source IPNS (Argonne Na- 
tional Laboratory, Argonne), and the cold-source triple- 
axis spectrometer SPINS at the NIST Center for Neutron 
Research. The time-of-flight measurements used incident 
energies of 6, 25, 35, and 60 meV, and were normalized 
on an absolute intensity scale using a vanadium standard. 
The higher energy runs showed no evidence of crystal 
field excitations above 20 meV. LRMECS has continu- 
ous detector coverage from 2.4° to 117.6°. The data at 
the highest angles (> 100°) are dominated by phonon 
scattering, but extrapolation of the momentum transfer 
dependence to low angle shows that the phonon contribu- 
tion below 30° is sufficiently small at all incident energies 
to be neglected in our analysis. The SPINS data, which 
were collected with a fixed final energy of 3.7 meV us- 
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ing a cold BeO filter, showed evidence of a well-resolved 
transition at 0.7 meV, which decreases in intensity with 
increasing temperature, i.e., it is a crystal field transition 
from the ground state. 

Figure ^ shows low-angle LRMECS data measured 
with an incident energy of 35 meV at 1.8, 3.2, 5, and 
20K. The data were summed from 2.4° to 30°, giving an 
average momentum transfer at the elastic position of 1.2 
A -1 . The data in Fig. Qshow that the magnetic spec- 
tra are dominated by two crystal field transitions; the 11 
meV peak is a ground-state transition, but the 17.2 meV 
peak only occurs at higher temperature and represents a 
transition from a low-lying excited state consistent with 
the 0.7 meV peak observed on SPINS. 



5 20 



T= 1.8K 




T = 5K 












A 




T = 3.2K 




T = 20K 








^^^^^^^ 


u I 


/ ^ \ 





-20 -10 10 20 -10 10 20 



Energy Transfer [meV] 

FIG. 1: S(Q,cj) measured at an average scattering angle of 
20° with an incident energy of 35 meV at four temperatures, 
normalized on an absolute scale. The average elastic momen- 
tum transfer is 1.2 A -1 . The solid line is the fit of the crystal 
field model described in the text with no adjustment of the 
intensity scale, and the dashed lines represent the individual 
crystal field transitions. 

The data are sufficiently complete to allow an unam- 
biguous determination of the crystal field potential. The 
praseodymium ions in PrOs4Sbi2 sit on lattice sites with 
cubic point group symmetry Th , which lacks two symme- 
try operations of the more common Oh group (C4 and 
C'2) UM- The crystal field levels in Th and Oh symmetry 
have the same degeneracies but different group theory 
labels and selection rules. In this article, we will use 
the more familiar cubic Oh labels, which can be mapped 
onto the Th labels using the following: Ti{Oh) — > T\(Th), 

r 3 (0 h ) -» r 23 (T h ), r 4 {o h ) -» r^(r h ), and r 5 {o h ) -» 

Th symmetry requires the addition of an extra crys- 
tal field parameter, B' 6 , to the conventional cubic crystal 
field Hamiltonian using the Steven's operator formalism: 



H CF = B 4 (O° 4 +5Oj)+B 6 (O 6 -21Ol)+B' 6 (O 2 6 -O 6 6 ) (1) 

The extra parameter, B' e , has a relatively small effect on 
the energies of the crystal field levels, but it mixes the 
T4 and T§ wavefunctions so that the usual cubic crystal 
field selection rules do not apply. Of particular impor- 
tance is that the dipole matrix element coupling the Ti 
singlet to the triplet is no longer zero. These dipole 
selection rules were a major factor in the previous crystal 
field assignments, and explain the discrepancy with the 
conclusions we outline below. 

In order to estimate of the importance of this term, we 
have performed a Superposition Model analysis of the B' e 
parameter for PrOs4Sbi2 0, 0] . The starting point of 
this model is the assumption that the crystal field is the 
superposition of two-body potentials due to the neigh- 
boring ligands. 

BY l = e l J2MRi)Kr(0i,<t> i ) (2) 

i 

where 0; are reduced matrix elements , Ai (R) are param- 
eters representing the strength of the two-body poten- 
tial, and K[""(d,(j)) are geometric functions tabulated in 
Ref. 0| . The sum is over all the neighboring ligands at 
(Ri,9i,4>i). The point-charge model obeys this superpo- 
sition principle but makes specific predictions concerning 
the values of Ai{R) that do not affect our conclusions. If 
we assume that the crystal field potential is dominated 
by the nearest-neighbor cage of twelve equidistant anti- 
mony ions, the model predicts the ratio of B' & to Bq with- 
out any adjustable parameters. Such a calculation shows 
that B'q is substantial in PrOs 4 Sbi 2 . B' % jB% = -53.4, so 
neglecting this term, as previously proposed |l2^. is not 
justified. 

In PrOs4Sbi2, the bulk susceptibility below 5K indi- 
cates that the crystal field ground state is non-magnetic. 
Therefore, the ground state is either the Y\ singlet or the 
T 3 doublet. Since the LRMECS data are normalized on 
an absolute scale and comprise measurements at a num- 
ber of temperatures, there are four ways in which the two 
models can be distinguished. 

a) Assuming that the 0.7 meV transition is the lowest- 
lying excitation, the absolute cross section of the 
T 3 — > T§ transition is considerably larger than the 
ri — > T$ transition and depends sensitively on the 
value of B' 6 . Although the latter transition is not 
dipole-forbidden in Th symmetry, it is still rela- 
tively weak. 

b) The absolute intensity of the other ground state 
transition at 11 meV is very different in the two 
models. The T 3 — > transition is 25% stronger 
than the Ti — ► T4 transition. 
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c) The intensities of the ground state transitions fall 
much more strongly with a Ti ground state than 
with a T 3 ground state because of the greater con- 
trast between the ground state and excited level 
degeneracies. For example, increasing the temper- 
ature from 1.8 K to 3 K is predicted to reduce the 
intensity of the 11 meV transition by 15% for the 
singlet ground state, but only 5% for the doublet 
ground state. 

d) There should be an excited state transition at 
~16.5 meV in the case of the Ti ground state. No 
such transition is predicted for a T3 ground state 
below 100 K. 

This gives us confidence that it is possible to estab- 
lish the crystal field level scheme unambiguously by using 
all the neutron data in a simultaneous refinement of the 
crystal field parameters, using the technique discussed in 
Ref. . We used the four spectra shown in Fig. and 
an additional spectrum at 10 K, which is not shown, with 
the additional constraint of requiring the lowest transi- 
tion to be at 0.7 meV. The peak lineshapes were Loren- 
zians convolved with the instrumental resolution. The 
linewidths did not vary significantly with temperature 
and were constrained to be equal in the final refinement. 
In the least-squares fitting procedure, the only adjustable 
parameters were the crystal field potential, the elastic in- 
tensity, and the common lincwidth. The reliability of the 
absolute normalization was such that we did not need to 
adjust the overall intensity scale. 
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FIG. 2: Temperature dependence of the intensity of the 11 
meV peak, represented as a Van Vleck susceptibility (see 
Ref. [Ig] for details), compared to the two possible crystal 
field models. 

We performed several refinements with starting param- 
eters consistent with both ground states. However, it was 
clear that only the Ti ground state was consistent with 
our data. As Fig. |3 shows, the temperature dependence 
of the ground state transitions is too strong to be repro- 
duced accurately with the T3 ground state. The intensity 



of the 0.7 meV peak is too weak in the 6 meV data (not 
shown). Furthermore, the T3 model does not reproduce 
the excited state transition at 17 meV observed above 
10 K. 

The results of the refinement are shown in Fig. ^ 
where all the transitions, from both ground and excited 
states, are plotted as dashed lines. It is clear that the Ti 
model is able to reproduce all the observed transitions. 
The crystal field potential has the following parameters: 
Bi = 0.20(1) x 10~ 2 meV, B 6 = 0.11(2) x 10~ 3 meV, 
and \B' 6 \ = 0.90(5) x 10~ 3 meV. The ratio, B' 6 /B 6 = 8.2, 
is smaller than predicted by the Superposition Model, 
which can be explained by adding a contribution to the 
crystal field potential from the osmium sublattice, which 
contributes to B 6 , but not to B' e . 

The conclusion that the Ti level is the crystal field 
ground state in PrOs4Sbi2 means that the quadrupo- 
lar Kondo effect cannot be responsible for the observed 
conduction electron mass enhancement. Nevertheless, it 
does not rule out other models in which the quadrupolar 
degrees of freedom of the rare earth f-electrons are im- 
portant. As we now discuss, the properties of PrOs4Sbi2 
can be explained by a delicate balance between two types 
of interaction, magnetic dipolar and quadrupolar, be- 
tween the conduction electrons and the praseodymium 
/-electrons. 

A theory of conduction electron mass enhancement due 
to inelastic scattering by crystal field transitions in a sin- 
glet ground state system was developed over twenty years 
ago jl^. According to Fulde and Jensen [l9|. the mass 
enhancement due to the inelastic transition at energy A 
between two levels, labeled i and j, is given by 

^ = l + (., 7 -l)^A(0)^M^ (3) 
mo J A 

where gj is the Lande factor, J s f is the exchange inte- 
gral coupling the conduction electrons to the f-electrons, 
and N(Q) is the bare conduction electron density-of- 
states at the Fermi level. If we assume that iV(0) is 
the same as in the isostructural lanthanum compound, 
then we can use the measured Sommerfeld coefficient 
in LaOs4Sbi2, 7 = 45 mJ/mol K 2 (averaging over the 
two published values 0, H3) in order to estimate that 
AT(0) = 37/27r 2 fc| = 9.6 x 10~ 3 meV _1 . We have 
no reliable estimate for J s f, but if we assume that the 
value derived in praseodymium metal, 0.085 eV [lflj . 
represents a reasonable order-of-magnitude, we obtain 
a mass enhancement of ~20. Given the uncertainty in 
the value of J s / , this is in reasonable agreement with ex- 
periment, falling between the estimates based on specific 
heat (~50 ||) and de Haas- van Alphen measurements 
(<7.6 0). 

The same s — f exchange that is responsible for the 
mass enhancement will tend to suppress superconductiv- 
ity through magnetic pair-breaking |2l| . However, the 
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superconducting transition temperature in PrOs4St>i2, 
T c = 1.85 K, is 2.5 times larger than the non-magnetic 
LaOs 4 Sbi2, T c = 0.74 K Y|. We propose that the reso- 
lution of this apparent discrepancy involves quadrupolar 
interactions that conserve time-reversal symmetry and 
therefore enhance pair formation ^5}. This effect, known 
as aspherical Coulomb scattering, is believed to reduce 
the rate of suppression of T c vs praseodymium concentra- 
tion in the singlet ground state system Lai-^Pr^Sna [52, 
[2^]. However, this would be the first compound in which 
T c is increased by praseodymium substitution [iij] . 

Fulde et al predict that aspherical Coulomb scatter- 
ing will produce the strongest enhancement of T c when 
A/T c is ~10 [l5j. In the case of Lai-^Pr^Sna, this ratio 
only occurs when T c has already been substantially sup- 
pressed by magnetic pair breaking |23| . However, if we 
assume that the crystal field potential is nearly constant 
in the Lai_ x Pr x Os4Sbi2 series, the optimum ratio occurs 
at x — > without requiring any suppression of T c |24j |. 
The strongest pair-breaking arises from the Ti — * T± 
transition, because it has the strongest dipole matrix ele- 
ments. This is at much higher energy (11 meV) than the 
Ti — > T$ transition (0.7 meV), which has a weak dipole 
but strong quadrupole matrix element, and is therefore 
responsible for the quadrupolar pair enhancement. In 
Lai_ ;r Pr a ;Sn3, these two transitions have comparable en- 
ergies. The crystal field level scheme in PrOs4Sbi2 is 
much more favorable for increasing T c through this mech- 
anism. 

In conclusion, we have performed comprehensive in- 
elastic neutron scattering measurements of the tem- 
perature dependence of the crystal field transitions in 
PrOs4Sbi2 which strongly suggest that the Ti singlet is 
the ground state. This would rule out the quadrupolar 
Kondo effect as the mechanism for the heavy fermion 
state, but favors another scenario in which the observed 
mass enhancement would arise from inelastic exchange 
scattering of the conduction electrons by the low-lying 
crystal field levels. We argue that inelastic quadrupolar 
scattering, also known as aspherical Coulomb scattering, 
provides an explanation for the enhancement in the su- 
perconducting transition temperature compared to the 
isostructural lanthanum compound. The importance of 
quadrupole interactions in PrOs4Sbi2 is evident in the 
antiferroquadrupolar order observed in high magnetic 
field [25J, |26| . Our results suggest that it plays a vital 
role in the superconducting phase as well. 
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